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Abstract 

In this paper we consider two models which exhibit equilibrium BEC superradiance. 
They are related to two different types of superradiant scattering observed in recent 
experiments. The hrst one corresponds to the amplification of matter-waves due to 
Raman superradiant scattering from a BE condensate, when the recoiled and the con- 
densed atoms are in different internal states. The main mechanism is stimulated Raman 
scattering in two-level atoms, which occurs in a superradiant way. Our second model 
is related to the superradiant Rayleigh scattering from a BE condensate. This again 
leads to a matter-waves amplification but now with the recoiled atoms in the same 
state as the atoms in the condensate. Here the recoiling atoms are able to interfere 
with the condensate at rest to form a matter- wave grating (interference fringes) which 
is observed experimentally. 
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1 Introduction 



This paper is the third in a series about models for equilibrium Bose-Einstein Condensation 
(BEC) superradiance motivated by the discovery of the Dicke superradiance and BEC matter 
waves amplification In these experiments the condensate is illuminated with a laser 

beam, the so called dressing beam. The BEC atoms then scatter photons from this beam 
and receive the corresponding recoil momentum producing coherent four-wave mixing of light 
and atoms The aim of our project is the construction of soluble statistical mechanical 
models for these phenomena. 

In the first paper [B], motivated by the principle of four-wave mixing of light and atoms [S], 
we considered two models with a linear interaction between Bose atoms and photons, one 
with a global gauge symmetry and another one in which this symmetry is broken. In both 
cases we provided a rigorous proof for the emergence of a cooperative effect between BEC 
and superradiance. We proved that there is equilibrium superradiance and also that there is 
an enhancement of condensation compared with that occurring in the case of the free Bose 
gas. 

In the second paper [Zj we formalized the ideas described in [H El by constructing a ther- 
modynamically stable model whose main ingredient is the two-level internal states of the 
Bose condensate atoms. We showed that our model is equivalent to a bosonized Dicke maser 
model. Besides determining its equilibrium states, we computed and analyzed the thermo- 
dynamic functions, again finding the existence of a cooperative effect between BEC and 
superradiance. Here the phase diagram turns out to be more complex due to the two-level 
atomic structure. 

In the present paper we study the effect of momentum recoil which was omitted in jUj and 
[7j. Here we consider two models motivated by two different types of superradiant scattering 
observed in recent experiments carried out by the MIT group, see e.g. PJ-jSj- Our first 
model (Model 1) corresponds to the Raman superradiant scattering from a cigar-shaped 
BE condensate considered in . This leads to the amplification of matter waves (recoiled 
atoms) in the situation when amplified and condensate atoms are in different internal states. 
The main mechanism is stimulated Raman scattering in two-level atoms, which occurs in a 
way similar to Dicke superradiance [Zj. 

Our second model (Model 2) is related to the superradiant Rayleigh scattering from a cigar- 
shaped BE condensate [2|, 0. This again leads to a matter- wave amplification but now with 
recoiled atoms in the same state as the condensate at rest. This is because the condensate 
is now illuminated by an off-resonant pump laser beam, so that for a long-pulse the atoms 
remain in their lower level states. In this case the (non- Dicke) superradiance is due to 
self-stimulated Bragg scattering [Sj. 

From a theoretical point of view both models are interesting as they describe homogeneous 
systems in which there is spontaneous breaking of translation invariance. In the case of the 
Rayleigh superradiance this means that the phase transition corresponding to BEC is at 
the same time also a transition into a matter-wave grating i.e. a "frozen" spatial density 
wave structure, see Section 4. The fact that recoiling atoms are able to interfere with the 
condensate at rest to form a matter-wave grating (interference fringes) has been recently 
observed experimentally, see (3]-(Ej, and discussion in [H] and (0]. 



In the case of the Raman superradiance there is an important difference: the internal atomic 
states for condensed and recoiled bosons are orthogonal. Therefore these bosons are different 
and consequently cannot interfere to produce a matter-wave grating as in the first case. 
Thus the observed spatial modulation is not in the atomic density of interfering recoiled 
and condensed bosons, but in the off-diagonal coherence and photon condensate producing 
a one-dimensional (corrugated) optical lattice, see discussion in Section 4. 

Now let us make the definition of our models more exact. Consider a system of identical 
bosons of mass m enclosed in a cube A C W of volume V = |A| centered at the origin. We 
impose periodic boundary conditions so that the momentum dual set is A* = {2ixp/V 1 ^ u \p E 
Z u }. 

In Model 1 the bosons have an internal structure which we model by considering them as 
two-level atoms, the two levels being denoted by a = ±. For momentum k and level a, a* k a 
and afc jCT are the usual boson creation and annihilation operators with [a*.^, a* k , a ,] = 5^8^. 
Let e{k) = \\k\\ 2 /2m be the single particle kinetic energy and N k>a = a* ka a k ^ the operator 
for the number of particles with momentum k and level a. Then the total kinetic energy is 

T 1A = ^e(£;)(iVV + AV) (1-1) 

and the total number operator is N 1>A = X^&eA*(-^M- + ^k,-)- We define the Hamiltonian 
H 1>A for Model 1 by 

H 1>A = T 1>A + U hA (1.2) 

where 

U hA = nb* q b q + -l=( a * q+ a _b q + a q+ a*_b* q ) + Ajv^, (1.3) 

g > and A > 0. Here b q , b* are the creation and annihilation operators of the photons, 
which we take as a one-mode boson field with [b q , b*] = 1 and a frequency Q. g is the coupling 
constant of the interaction of the bosons with the photon external field which, without loss 
of generality, we can take to be positive as we can always incorporate the sign of g into 
b. Finally the A-term is added in (jl.2)l to obtain a thermodynamical stable system and to 
ensure the right thermodynamic behaviour. This is explained in Section 2. 

In Model 2 we consider the situation when the excited atoms have already irradiated photons, 
i.e. we deal only with de-excited atoms a — — . In other words, we neglect the atom excitation 
and consider only elastic atom-photon scattering. This is close to the experimental situation 
in which the atoms in the BE condensate are irradiated by off -resonance laser beam. 
Assuming that detuning between the optical fields and the atomic two-level resonance is 
much larger that the natural line of the atomic transition (superradiant Rayleigh regime 
[21 E]) we get that the atoms always remain in their lower internal energy state. We can 
then ignore the internal structure of the atoms and let a* k and a k be the usual boson creation 
and annihilation operators for momentum k with [afe,a^/] = 5k t k', Nk = a* k ak the operator 
for the number of particles with momentum k, 

T 2 , A = < k ) N * ( L4 ) 

fceA* 

the total kinetic energy, and A^ 2jA = YlkeA* ^k the total number operator. We then define 

(2) 

the Hamiltonian H K for Model 2 by 

H 2tA = T 2A + U 2iA (1.5) 



where 

U 2 , A = b* q b q + ^Ka 6 g + a ? a£6*) + ^iV 2 2 A • (1.6) 
This paper is structured as follows: 

In Section 2 we give a complete rigorous solution of the variational principle for the equi- 
librium state for Model 1 and we compute the corresponding pressure as a function of the 
temperature and the chemical potential. We prove that this model exhibits Raman super- 
radiance. 

In Section 3 we study Model 2 and show that Rayleigh superradiance occurs in this case. 
The analysis is very similar to that of Model 1 and therefore we do not repeat it but simply 
state the results. 

In Section 4 we show that in both models there is spontaneous breaking of translation invari- 
ance in the equilibrium state. We relate this with the spatial modulation of matter-waves 
[matter-wave grating). We find that in Model 1 there is no such spatial modulation in spite 
of the breaking of translation invariance while in Model 2 this spatial modulation exists. We 
conclude with several remarks. 

We close this introduction with the following comments: 

- In our models we do not use for effective photon-boson interaction the four-wave mixing 
principle, see jS], jEj, jlt)j . The latter seems to be important for the geometry, when a lin- 
early polarized pump laser beam is incident in a direction perpendicular to the long axis 
of a cigar-shaped BE condensate, inducing the "45°- recoil pattern" picture 0-0. Instead 
as in [7], we consider a minimal photon-atom interaction only with superradiated photons, 
cf [U]. This corresponds to superradiance in a "one-dimensional" geometry, when a pump 
laser beam is collimated and aligned along the long axis of a cigar-shaped BE condensate, 
see 0, [12]. 

- In this geometry the superradiant photons and recoiled matter-waves propagate in the same 
direction as the incident pump laser beam. If one considers it as a classical "source" (see 
0), then we get a minimal photon-atom interaction [Zj generalized to take into account the 
effects of recoil. Notice that the further approximation of the BEC operators by c-numbers 
leads to a bilinear photon-atom interaction studied in pj, (Sj- 

- In this paper we study equilibrium BEC superradiance while the experimental situation (as 
is the case with Dicke superradiance [T^]) is more accurately described by non-equilibrium 
statistical mechanics. However we believe that for the purpose of understanding the quantum 
coherence interaction between light and the BE condensate our analysis is as instructive and 
is in the same spirit as the rigorous study of the Dicke model in thermodynamic equilibrium, 
see e.g. |E|-P5|. 

- In spite of the simplicity of our exactly soluble Models 1 and 2 they are able to demon- 
strate the main features of the BEC superradiance with recoil: the photon-boson condensate 
enhancement with formation of the light corrugated optical lattice and the matter-wave grat- 
ing. The corresponding phase diagrams are very similar to those in However, though 
the type of behaviour is similar, this is now partially due to the momentum recoil and not 
entirely to the internal atomic level structure. 



2 Model 1 



2.1 The effective Hamiltonian 



We start with the stability of Hamiltonian (J1.2|) . Consider the term ?7 liA in ()1.3j) . This gives 



^-^S + aoJ(^ + ^w^ a >o-) - ^yN .(N q+ + 1) + — „ ljA 
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£ 2^ iV "-4lW iV °-' iV - + 1 »- (") 

On the basis of the trivial inequality Aab < (a + b) 2 , the last term in the lower bound in (j2.1j) 
is dominated by the first term if A > g 2 /8Q, that is if the stabilizing coupling A is large with 
respect to the coupling constant g or if the external frequency is large enough. We therefore 
assume the stability condition: A > g 2 /8Q. 

Since we want to study the equilibrium properties of the model (jl.2)l in the grand-canonical 
ensemble, we shall work with the Hamiltonian 

H liA (ji) = H liA -fiN ltA (2.2) 

where \i is the chemical potential. Since T 1A and the interaction U ljA conserve the quasi- 
momentum, Hamiltonian (jl.2j) describes a homogeneous (translation invariant) system. To 
see this explicitly, notice that the external laser field possesses a natural quasi-local structure 
as the Fourier transform of the field operator b(x): 

b q = [ dxe iq - x b(x). (2.3) 



A 



If for z G W, we let r x be the translation automorphism (r z b)(x) = b(x + z), then since 
we have periodic boundary conditions, r z {b q ) = e~ iqz b q and similarly r z (afc j(J ) = e~ %kz a^^. 
Therefore, the Hamiltonian (jl.2)l is translation invariant. Consequently, in the thermody- 
namic limit, it is natural to look for translation invariant or homogeneous equilibrium states 
at all inverse temperatures f3 and all values of the chemical potential ji. 

Because the interaction (J1.3)) is not bilinear or quadratic in the creation and annihilation 
operators the system cannot be diagonalized by a standard symplectic or Bogoliubov trans- 
formation. Therefore at the first glance one is led to conclude that the model is not soluble. 
However on closer inspection one notices that all the interaction terms contain space aver- 
ages, namely, either 

\- [ dxa_(x), (2.4) 



W v JA 

and its adjoint, or 

— / dxa* a (x)a (T (x). (2.5) 



A 



Without going into all the mathematical details it is well-known that space averages 
tend weakly to a multiple of the identity operator for all space-homogeneous extremal or 
mixing states. Moreover as all the methods of characterizing the equilibrium states (e.g. the 
variational principle, the KMS-condition, the characterization by correlation inequalities etc. 



[17]) involve only affine functionals on the states, we can limit ourselves to looking for the 
extremal or mixing equilibrium states and in so doing we can exploit the above mentioned 
property for space averages. One way of accomplishing this is by applying the so-called 
effective Hamiltonian method, which is based on the fact that an equilibrium state is not 
determined by the Hamiltonian but by its Liouvillian. The best route to prove the exactness 
of the effective Hamiltonian method (cf ^Hj) is to use the characterization of the equilibrium 
state by means of the correlation inequalities ^H], [T7j : 

A state uj is an equilibrium state for H 1 . A (/i) at inverse temperature (3, if and only if for all 
local observables A, it satisfies 

lim 0u> ([A*, [H 1>A (p), A]}) > u(A*A) In (2.6) 

Clearly only the Liouvillian [H lyA (p), ■} of the Hamiltonian enters into these inequalities and 
therefore we can replace H 1A (p) by a simpler Hamiltonian, the effective Hamiltonian, which 
gives in the limiting state u the same Liouvillian as H 1A (p) and then look for the equilibrium 
states corresponding to it. Now in our case for an extremal or mixing state uj we define the 
effective translation invariant Hamiltonian H° s A (p,r], p) such that for all local observables A 
and B 

lim u (A, [H h M, B}) = lim cu (A, [H*(ji, 77, p), B}) . (2.7) 
The significance of the parameters r] and p will become clear below. One can then replace 

by 

Km Pu,([A\[H*fa V ,p),A]]) > w (AM)ln|^^. ( 2 - 8 ) 

We can choose H^ a A (p,r], p) so that it can be diagonalized and thus ()2.8|) can be solved 
explicitly. For a given chemical potential p, the inequalities ()2.8|) can have more than one 
solution. We determine the physical solution by minimizing the free energy density with 
respect to the set of states or equivalently by maximizing the grand canonical pressure on 
this set. 

Let the effective Hamiltonian be defined by 

H iA^V,p) = (Xp - p + e(q))a* + a q+ + {Xp - p)a* _a - + ^(va* q+ b q + fja q+ b*) 

+n b* q b q + ^ (Coo- + (a* _) + K A + (Ap - p)N' hA (2.9) 

where 

T U= E E ( 2 - lQ ) 

N U= E E N *,-> ( 2 - n ) 

r) and ( are complex numbers and p is a positive real number. Notice that the Hamiltonian 
()2.9|) is translation invariant, but it is not gauge invariant for ( 7^ because of the linear 
terms in aQ_,ao_ operators. Therefore in this case, ()2.9|) generates translation invariant 
states, which are extremal with respect to the gauge group in the zero-minus mode. They 
are labeled by the arg (. One can easily check that (|2.7j) is satisfied if 

V= -j= . C = y and p = — - — , (2.12) 



where the state ui coincides with the equilibrium state ( ■ )Hf (n,v,p) defined by the effec- 
tive Hamiltonian H° a A (fi,r], p). From (|2.7|) and (J2.12J) we then obtain the self- consistency 
equations 

1 1 1 

(a -) ff 8ff ( = —(aq+b*)^^^ ^, p = — (iV ljA ) H eff ^ p ). (2.13) 



y \-u-/ ni~ A yi,ri,p)-> •> y \-y-r-g/ n^ A (fi,ri,p)i r y 

The structure of ()2.9j) implies that the parameter ( is a function of p and p through ()2.13j) . 
So, we do not need to label the effective Hamiltonian by The important simplification 
here is that H° s A (p,r], p) can be diagonalized: 

H i G a(PiViP) = E + (p,r],p)a* 1 a 1 + E^(p,r],p)a* 2 a 2 (2.14) 

a 2 V\C\ 2 

+(\ P - p) a ;a 3 + r 1A + (x P - p)N' hA + * , 

where 

1 1 

E + (p,r],p) = ~(n-p + \p + e{q)) + - v 7 ^ + P - Ap - e(g)) 2 + # 2 |p| 2 , 

E_(p,r),p) = l -(Sl-p + \p + e(q)) -^(tt + p-Xp- e(g)) 2 + ^F, 

(2.15) 



«i = a g+ cos# + b q sin 9, a 2 = a q+ sin 9 — b q cos 9, a 3 = a - + 7777 r? (2.16) 

2(Ap - p) 

and 

tan20 = -- ^ --. (2.17) 

tt + p-Xp- e(q) 

Note that the correlation inequalities (J2.8|) (see [TH]) imply that 

lim uj (A*, [H 1A (p), A]) > (2.18) 

V^oo 

for all observables A. Applying (J2.18)) with A = a^,, one gets the condition Xp — p > 0. 
Similarly, one obtains the condition Ap + e(g) — p > by applying (|2.18j) to A = a* + . We also 
have that E + (p,r),p) > EL(p, p,p) and E_(p,rj,p) = when |p| 2 = 4f2(Ap + e(q) — p)/g 2 
and then E + (p, p, p) = Q — p + Ap + e(q). Thus we have the constraint: 

\r 1 \ 2 <AQ{Xp + e{q)-p)/g 2 (2.19) 

We shall need the above information to make sense of the thermodynamic functions below. 
Of course the parameters p, p and consequently E± are V dependent but for simplicity we 
do not indicate this dependence explicitly. 

We can foresee that for some values of p there will be Bose-Einstein condensation in the mode 
{0—}. We know that in this case the gauge invariant, homogeneous states are not extremal 
within the class of translation invariant equilibrium states [20] • Therefore to ensure that the 
states that we shall obtain are extremal we add to the Hamiltonian a gauge breaking term 



(ha _ + ha* _) , (2.20) 



and then let h — > after the thermodynamic limit. The corresponding effective Hamiltonian 
is then 

H*Qm, r), p\ h) := H*(jm, V,P)~^ (ha - + K-) (2-21) 
with h G C. The equations corresponding to (J2.13|) now become 

^ = ^/^( a 0-)ff=ff (ji,r),p;h)i C = y ( a q+Ki)H°f A (p,r), P ;h)i P = y{Ni, A ) Hff A (/j,,ri, P ;h)- (2.22) 

These consistency equations can be made explicit by using the above diagonalization: 

^W^P) {C - kl (2 - 23) 

C= I ill I I \ \ (2 24) 

* 2 V(E + - EJ) \ eP E + - 1 eP E - - 1 J 1 ' 



and 



,2 1 1 1 



^ + V e-^~ x P) - 1 + 2V I eP E + - 1 ^ 
(fi-Xp- e(q) + ft) 



2V r (£+ - £L) [ - 1 e? E - - 
Combining ()2.23j) and (J2.24|) we obtain the equation: 

= sfv ( i i 1 gfr (226) 

V 4(p - Xp)V(E+ - E.) \ e? E + - 1 eP E - - 1 J 2(ji-\p)' { ' ' 

It is now clear that the equilibrium states are determined by the limiting form of the con- 
sistency equations (|2.23j) - (|2.26j) . We solve these equations and obtain the corresponding 
pressure so that we can determine the equilibrium state, when there are several solutions for 
a particular chemical potential. 

We shall need the following definitions: 

= Wf L/ k J^T- (2 ' 27) 

and 

poQi) = f d 3 fcln(l-e-^ e(fc) -^), (2.29) 



(2vr)3 , 

that is the grand-canonical energy density, the particle density and the pressure for the free 
Bos e- gas for p < 0. Let 

s (p)=(3(e (p)+p (p)), (2.30) 

and note that Sq(p) is an increasing function of p. We shall denote the free Bose-gas critical 
density by p c , i.e. p c := po(0). Recall that p c is infinite for v < 3 and finite for v > 3. 



2.2 Solution of consistency equations 



Notice first that equations (J2.23)) and (J2.24j) imply 

arg ( = arg r] = arg h = <p , (2-31) 

i.e. one can consider the corresponding parameters in the consistency equations (|2.23J) - 
()2.26j) to be real and non-negative. 

Remark 2.1 By virtue of the upper bound h2.19i) and equation 112.26]) we get that 5 := 
Xp — p > for any volume V as soon as h ^ 0. By the same reasoning one gets from 
\2. 2b}) that Ymvy^^ inf E_ > ( Case A ), or at most limy^oo EL = in such a way that 
limy^oo VE_ is finite (Case B). 

We start our analysis of the solution of the consistency equations from small densities (small 
chemical potentials), when there is no condensates, passing then to higher values. So, later 
on we distinguish the two possibilities indicated in Remark 12. II 

Case A : lim > 0. (2.32) 

V— >oo 

Then by Remark 12. II the consistency equations (J2.23j) - (J2.26)) in the thermodynamic limit 
yield: 

V = ^ , C = 0, (2.33) 
and the equation for the particle density (J2.25)) takes the form: 

p = 2\ Po (-5)-5+^f. (2.34) 

We also have the limiting expectations 

1 2 g 2 \h\ 2 

y 1 ™ ^(W^W^N = (2 ' 35) 

and (by virtue of (j232)) 

y 1 ™ ^K+Vf>Hf A(/ v,,M) = ^ y(K b «)H°e A ( m h) = 0- (2.36) 

Let us now examine equation ()2.34j) . For h ^ it is solvable for any p and we shall denote 
its unique solution by S(p, h) = Xp(p, h) — p. 

Solution 1: Suppose that p < p c := 2\p c = 2Apo(0). Since \h\/5 —> as h — > 0, then 
5(p, h) — > S(p), where 5(p) is the unique solution of equation: 

p = 2Xp (-5) - 5. (2.37) 

Then we see from ()2.35|) and ()2.36|) that in this case in the thermodynamic limit there is no 
condensation in the {0—} and other modes: 

J™ y K- a o->fff A ( M ,o,P) = y L ™ = ^Wnf^M = 0- (2.38) 



After letting h — > 0, the energy density is given by 

lim — 

v^oo V 



lim V (H h M) HfA{llAp) = 2e {-S{n)) - 2{6{n) + n)po{-6{n)) + ~A(5(m) + fif 



= 2e M )-\m±J^ (2 .39) 
and the entropy density is equal to 

s(/x) = 2s (-5(/i)). (2.40) 
Since the grand-canonical pressure is given by 

then 

p(Ai) = 2p (-5(/i)) + ^ (W A +/i)2 . (2.42) 

Solution 2: Now suppose that fi > fi c = 2Xp c . Then to verify equation (J2.34)) in the limit 
h — > the solution must converge to zero: 5(/i, h) — > 0, in such a way that 

\a 2 \h\ 2 

w^k " - 2Aft - (2 - 43) 

Therefore it follows from (J2.34)) and (|2.H5jl that 
and (again by (|2.32|l ) the limit ()2.36|) gives 

i.e., there is no condensation in the q ^ modes and the laser boson field. 
In this case the energy density is given by: 

& ^<^M>ii- (*,,,,) = 2 ^o(0) - ^ (2.46) 
and the entropy density has the form: 

s(ji) = 2s (0) = 2(3(e (0) + Po (0)). (2.47) 
Thus for the pressure one gets: 

p(ji) = 2p (0) + ^ (2.48) 

Notice that the bound (|2.19j) and 1)2.44)1 imply the upper limit on chemical potential 

4^Ae(g) , A , 

/i < /i c + ^ 2.49 



for which Solution 2 applies. 

This means that for the higher densities or chemical potentials: 

/u>/i c + r^, 2.50 

a 2 

to satisfy equation (|2.25|) we have to consider 

Case B : lim = 0. (2.51) 

V— >oo 

Then (see f!2.19|) and Remark 12. lj) in the thermodynamic limit: |?7| 2 — > 4Q(5 + e(q))/g 2 . In 
fact, to obtain a finite limit in (J2.26)) the corresponding large- volume asymptotic should to 
be 

N'-F^^-^j (2 ' 52) 

for some r > 0. This implies that 

E^n + <q) - S , E.» ? (2.53) 



and (J2. 2(1)1 becomes in the limit: 



v [1-£l)=!±. (2.54) 
45fi/ 26 v ; 



The last equation gives 

T = ^_M. (2.55) 
9V 

Taking the limit V — > oo in ()2.25|) we get 

4Afi(5 + e(g)) 4A5fi 2A/iO nx . . 

= - + — + 2Xp (-5) - 5. 2.56 

9 2 9 9V 

We can also check by using the diagonalization that in this case 

,- 1 / . v I i2 4^ + 6(g)) 

1 / , Asn 2hn . 

^v^ aq+aq+ ' H ^ {mh) = T = 1 2 ~~ (2 ' 58) 

Bmiw^^/i) = ^ = + (p " ■ (2 ' 59) 

We note here that if we take 5 = in (J2.56)) . then for h — > this expression gives the limiting 
value of the chemical potential ()2.49|) for Solution 2. 

Solution 3: Let fi > fj, c + 4f2Ae(g) / 'g 2 , see 1)2 .50|) . Now we take Case B and let h — > 0. 
Then by 1)2.57)1 . (J2.58)) and (J2.59)) we obtain a simultaneous condensation of the excited/non- 
excited bosons and the laser photons in the g-mode: 

1 / * l2 4fi(5 + e(g)) 
^ OD y( a o- a o-)H^,r,, P ) = \v\ = ~ 2 , (2-60) 



Mm ± W*-^, = (2.62) 



Equation (|2.5fj|) becomes: 



4Xn(5 + e(q)) 4\sn ^ 

fJ>= 5 1 5— + 2Xp (-b) - 6. (2.63) 

9 9 



Using the diagonalization of (J2.9J) one computes also 



vfe yK+^ftOw) = ~2 • ( 2 - 64 ) 

In this case the energy density is given by: 

1 A5il An(5 + e(q)) 805(5 + e(q)) 

hm — (# 1)A (/i))tfeff A(/w) = (e(g) -/i)— ^- -/i- 



2 



4(5 + e(g))<5 _ . . . 1 (5 + p) 2 
+^ ^ 2 + 2£ o(-5) - 2(5 + /i)poM) + - V A P; 

= "a + 2£ o{-o)-- ^ • (2.65) 

The entropy density is again given by 

s(p) = 2s (p - Xp) (2.66) 

and the pressure becomes 

l(5 + p) 2 4Q(5 + e(q))5 
A 9 2 



PW = 2p (-5) + . (2.67) 



Notice that only \rj\ is determined and not the phase (J2.31|) of rj , that is that we can only 
determine the state up to a gauge transformation in the a = — fields. 

The fact that here we have condensation in a state with non-zero momentum is extremely 
significant and is related to the spontaneous breaking of translation invariance. We shall 
examine this important aspect in Section 4. 

The Solutions 1, 2 and 3 represent possible equilibrium states. For a given value of the 
chemical potential p, two or even three of these solutions may be possible, see Figure 1 
below. To distinguish between them we have to compare the corresponding pressures to 
determine which is maximum. The analysis, which is given in the next subsection, involves 
a detailed study of the pressure. We find that the situation is as described below. 

Let k = 8QX/g 2 — 1 and a = e(q)(n + l)/2. From the condition for thermodynamic stability 
we know that k > 0. In this notation Solution 2 applies for p c < p < p c + a. Let 5q be 
the unique value of 5 G [0, 00) such that 2Xp' (— 5) = k and let p = 2Xp (—5 ) + k5 . Note 
that p < 2Xp c . 

The case when p + a > 2Xp c is easy. In this situation Solution 1 applies for p < 2Xp c 



and there exists pi > po + a (see definition after ()2.8()j) ) such that Solution 2 applies for 
2Ap c < p < \i\ and Solution 3 for p > 

When po + o < 2Ap c the situation is more subtle. In Subsection 2.3 we shall show that there 
exists fii > fio + a (I2.8()jl . such that Solution 3 applies for p > pi. However we are not 
able to decide on which side of 2Ap c , the point pi lies. If pi > 2Ap c the situation is as in 
the previous subcase, while if p + ol < pi < 2Ap c the intermediate phase where Solution 2 
obtains is eliminated. This the situation is similar to [7j, where one has a = 0. 

Note that for v < 3, Solution 1 applies when p < \i\ and Solution 3 when p > pi. 



2.3 The Pressure for Model 1 



This subsection is devoted to a detailed study of the pressure for Model 1 as a function of 
the chemical potential p. 

Recall that 5 is the limiting value of Ap — p, k = 8QX/g 2 — 1 and a = e(q)(n + l)/2. From 
above we have the following classification: 

Solution 1: Here p < p c . The density equation 

p = 2Xp (-5) - 5 (2.68) 

has a unique solution in 5, denoted by 5i(p) (previously denoted by 5(p)). Let 

p 1 (5,p):=2po(-5) + ^^. (2.69) 

Then 

p(ji)=Pi(S 1 (ji),fi). (2.70) 
Solution 2: Here p > p c , 5 = and the pressure is given by 

p(p)=p 2 (p):=2p (0) + ^. (2.71) 

Solution 3: The equation (|2.63|) can be re-written as 

p = 2\p Q (-5) + k5 + a. (2.72) 

Recall that S is the unique value of 5 G [0, oo) such that 2Xp' (— 5) = k, p = 2Ap (— S ) + kS 
and that p < 2Ap c . 

Then for p < p + a, equation (J2.72|) has no solutions. For p + a < p < 2Ap c + a 
this equation has two solutions: 5 3 (p) and <5 3 (p), where (5 3 (p) < 5 3 (p) if p 7^ p + a, and 
^3(^0 + a ) — ^3(^0 + «)• Finally for p > 2Ap c + a it has a unique solution <5 3 (p). Let 

„(,, „, : . 2po( _,) + W + -(«+!)* ,2.73) 



Then 



dp 3 (S 3 (p),p) = S 3 (p) + p < £ 3 (p) +P = dp 3 (5 3 (p),p) 
<ip A A <ip 



for p ^ p + a. Since p 3 (5 3 (po + a), Po + ") = P3(5 3 (p + a), Po + a), 

MW,IJ>)<P3(6 3 (iM),it) (2.75) 
for p + a < p < 2\p c + a. Therefore 

p(^)=P3(S 3 (ji),ti) (2.76) 

for all p > Hq + a. 

Note that S 3 (2Xp c + a) = so that 

p 3 (5 3 (2Ap c + a), 2Ap c + a)= p 1 (5 l (2Xp c ), 2\ Pc ) = p 2 (2Ap c ) = 2p (0) + 2\p\. (2.77) 
therefore 

^(^(po + a), /i + «) = ^(^(po + a), /i + «) < 2p (0) + 2Ap^. (2.78) 

Also for large p, p 3 (<5 3 (p),p) « (p 2 /2A)((k + l)/«) while P2(p) ~ (/ i2 /2A), so that 
^3(^3 (A*); A*) > P2(p) eventually. We remark finally that the slope of p 3 (<5 3 (p), p) is greater 
than that of P2(p) ; 

dp3{S 3 { P ),fi) = 5M + > /£ = ^(p) , 2 7g ^ 

dp A A dp ' 

so that the corresponding curves intersect at most once. 
The case a = 0, i.e. e(g = 0) = 0, has been examined in jjj. 
For the case a > we have two subcases: 

The subcase po + ol > 2Ap c is easy. In this situation Solution 1 applies for p < 2Ap c . From 
we see that 

P3(5 3 (A*o + «)> A*o + ol) < 2p (0) + 2Xp 2 c < p 2 (p + a (2.80) 

and therefore from the behaviour for large p we can deduce that there exists pi > po + ol 
such that Solution 2 applies for 2Ap c < p < p\ and Solution 3 for p > fi 1 . 

The subcase po + a < 2Ap c is more complicated. In FigureHJwe have drawn y = 2Xp (—5) — 5 
and y = 2Apo(— 5) + k5 + a for this subcase. We know that 

p 3 (5 3 (2Ap c ), 2Ap c ) < p 3 (5 3 (2Ap c + a), 2Ap c + a) = p 1 {S 1 {2Xp c ), 2Xp c ). (2.81) 

Therefore since the slope of p 3 {5 3 {pL), p) is greater than the slope of pi(<5i(p), p) for p + a < 
p < 2Ap c , (see Figure HJ): 

dps^M^ = S 3 {p) + p > ^i(p) + /i = rfp!(^(p),p) g2 
dp A A dp 

we can conclude that 

P-i(S 3 (fi + a),p + a) = p 3 (<53(Ak) + a),Po + a) < Pi(5i(p + a),Po + a). (2.83) 

We also know by the arguments above that there exists pi > po + ol such that Solution 
3 applies for for p > pi. However we do know on which side of 2Ap c , the point pi lies. 
If pi > 2Ap c the situation is as in the previous subcase while if p + a < pi < 2Ap c the 
intermediate phase where Solution 2 obtains is eliminated. 




Figure 1: The density equation for h = and v > 3 

3 Model 2 



As we said in the introduction the analysis for this model is very similar to that of Model 1. 
Therefore we briefly summarize the results without repeating the details. For Model 2 the 
effective Hamiltonian is 

#2,aO ; V, p) = (Ap - p + e(q))a* q a q + (Ap - p)a*a + |(p<& 9 + fja q b* q ) 

+0 b* q b q + ^ (Ca + (a* ) + T 2 ' iA + (Ap - p)^ (3.1) 

where 

The parameters p, ( and p satisfy the self-consistency equations: 

P = ^<a } H e ffA(/w) , C = y (a* q b q )Hf A („, V>P ), P = y( N **)H$Sr,r,, P y ( 3 - 4 ) 



Using the external sources (|2.2()jl and the same treatment as for the Model 1 in Section 2.2, 
we again obtain three cases: 

Solution 1: /i < Xp c = \i c . In this case there is no condensation: 

y 1 ™ y( a O a v)H^, V , P ) = ^yi^l^to) = ^ yWH^i^p) = 0, (3-5) 

the density equation is 

H = \po(S) - 6 (3.6) 

and the pressure is 

p{p)=Po{-5) + \^±^. (3.7) 



Solution 2: Let /i c < fi < fi c + 4ftAe(g) / g 2 . Then 5(p:, h) = lim /l ^ (^p(/ i ) h) — p) = 0. 

I 7 '' 2 = J im 77( a o a o)Hf.(v, v , P ) = T - Pc- (3i 



There is condensation in the k = mode but there is no condensation in the k = q mode 
and of the photon laser field: 

feyWOw) = y <W*!Ta0m7,p) = °- ( 3 - 9 ) 

The pressure density is given by 

p(Ai)= Po(0) + ^. (3.10) 

Solution 3: Let \x > \x c + 4ftAe(g) / g 2 , see (|2.50|) . Then there is simultaneous condensation 
of the zero-mode and the g-mode bosons as well as the laser g-mode photons: 



4ft (5 + e(q)) 
9 2 



1 . „ . 45ft . 

y 1 ™ y^WH^^p) = —T, (3.12) 

r 1 //,*/, \ 4(£ + e(g))£ 

J™ 7 <W^>,^) = ^ > ( 3 - 13 ) 



The density equation is 



r 4Vftg+gg) 



8 Aft 

l M=—(5 + e{q)/2)+p (-S)-8 (3.15) 



and pressure is 

l (5 + p) 2 4Q{5 + e{q))5 
2 A g 



P(W =Po(-0) + o x la • I 3 " 16 ) 



Note that relations between the values of p, in the three cases above are exactly the same as 
for Model 1 apart from the fact that 2p is now replaced by p an d 2p c by p c . To see this 
one has to compare the kinetic energy operators (jl.lj) and (|1.4j) . 



4 Spontaneous Breaking of Translation Invariance and Matter- Wave Grating 



The recently observed phenomenon of periodic spatial variation in the boson-density is 
responsible for the light and matter-wave amplification in superradiant condensation, see |2j- 
[3], ^21- This so called matter-wave grating is produced by the interference of two different 
macroscopically occupied momentum states: the first corresponds to a macroscopic number 
of recoiled bosons and the second to residual BE condensate at rest. Clearly this cannot 
happen in the translation invariant states and so it must be due to a spontaneous breaking 
of this invariance. 

Let us consider the situation in Solution 3 for Model 1. For simplicity we shall take q = 
(27r/7)e!, with e! = (1, 0, . . . , 0) G W and 7 > and we shall denote the limit Gibbs state 
for the effective Hamiltonian by uj\ 

We know that the existence of condensation in the zero-mode of the a = — bosons implies 
that the extremal states have broken gauge symmetry in the corresponding fields in this 
mode. As was remarked earlier this is indicated here by the fact that 77 is not zero. 

It has been shown in [20] that condensation in the zero-mode implies that the gauge invariant 
spatially homogeneous equilibrium states are not extremal but can be decomposed into a 
convex combination of extremal space homogeneous equilibrium states with broken gauge 
symmetry (spontaneous braking of gauge symmetry). 

From (J2JHJ), flZEZD and (gUm we see that in Solution 3 



lim 

y-»oo 



lim — u(a* + a q+ ) lim -u;(b*b q ). (4.2) 



This strongly suggests a similar decomposition when there is condensation in the g-mode. 
In fact when there is condensation in a mode g ^ one can argue, along the same lines 
as in [20], that the spatially homogeneous equilibrium states when considered within the 
equilibrium states which are periodic in the ei direction with period 7, are not extremal. 
They can be decomposed into a convex combination of extremal periodic equilibrium states 
which are not spatially homogeneous: 

1 P 

u = — dx u) x (4.3) 
7 Jo 

where 

Wx r 7ei = uj x (4.4) 

and 

(4.5) 

Therefore in this model we have spontaneous breaking of translation invariance. The rigorous 
and explicit construction of the states u x involves mathematical details which are outside of 
the scope of the present paper and is carried out in 



Let Ai = {x\x G A, < x < 7} and let V\ = |Ai|. Then we can write, for example, 

1 N 1 r 

lim a*/W= lim — V — / e ±iqx a*(x)dx. (4.6) 



Therefore in the representation corresponding to each of the extremal states u x , af + /VV 
converges weakly to a complex number which by (|4.4jl and ()4.5|) is equal to 



e Tiqx lim 



u {a 



* ) 

q+)i 



V^oo y'y 

where u = lo x=0 . It then follows from (|2.61j) . with 5 = limy^ QO (Ap — p), that 



(4.7) 



1 



,J il ^_77 a; °K+ a 9+) 



1 



lim 



45Q 



(4.1 



Similarly af_/y/V and b#/y/V converge weakly to a complex numbers and from (J2.60P and 
(I2.62j) respectively we obtain 



and 



lim — uj (a* _a „ 

V— >oo V 



Jim ^juj (b*b q ) 



lim —= 

v^oo Jy 



V^oo V 



The weak convergence of af + /VV and Oq_/\/V to complex numbers also implies that 



lim —= 

,# 



u (b q 



4Q(5 + e{q)) 



(4.9) 



(4.10) 



lim -u (a* + a k , + ) = lim -^=u {a* + ) lim -^=u {a k , + ) = 

V— +oo V V— >oo yy V— too yy 



fork^q and 



lim 2-uo(aZ_a k) -) = lim -^=u (a* ) lim -^=u (a kr 

V— >oo 1/ y— »oo yy V-+OC yy 



(4.11) 



(4.12) 



for k + 0. 

An alternative strategy to the one we have developed above would be to use the traditional 
method of introducing source terms in the Hamiltonian (jl.2j) to break both the gauge and 
translation symmetries and let: 

H 1A (0 := H 1>A - W(Zb* q + £b q ) . (4.13) 

Then the effective Hamiltonian becomes 

#T,a(/^ P\ C, = (A/0 - H + e(q))a*a q+ + (Xp - p)a* _a ^ 
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+^ i (C + 7p a >o- + (C + hK-°<7+ j- + n 6*6. 



l^l 2 1 

+^|C| 2 - ^ + T[ A + (Ap - /i)^, - -Vp 2 , (4.14) 



where 



b q = b q ->/V(C + y 



(4.15) 



and we can carry out the procedure of Section 2.1 to obtain the results equivalent to (|4.8|) - 



We now want to examine the possibility of interference between two different macroscopically 
occupied momentum states in Model 1 in the periodic states. Without loss of generality we 
can restrict ourselves to the state Uq. In this state the mean local particle density for the 
a = + bosons is 



p+(x) = lim . 



(k—p)x 



uj (a k+ a P!+J 



p + lim — 



fceA* peA* fceA* peA* p^fc 

(4.16) 

We know that condensation occurs only in the g-mode for the a = + bosons and in the 
q — - mode for the a = — bosons and therefore only the terms containing q survive in the 
integral sum ()4.16|) in the thermodynamic limit: 



P+K x ) 



p + lim - V 2 { 



e« k -^u (a* k>+ a q , + )} 



(4.17) 



fceA* k^q 



From the above discussion, in particular from (|4.11jh we see that p+(x) = p. Similarly 
P- (x) = p. Thus in spite of the fact that the state u is not space homogeneous, the total 
particle density is constant and equal to 2p. This means that in Model 1 we get no particle 
density space variation even in the presence of the light corrugated lattice of condensed 
photons, see (14.101) . 



Let us now turn our attention to the corresponding situation for Model 2 in Solution 3. The 
decomposition into periodic states still stands and again we have spontaneous breaking of 
gauge symmetry. In the representation corresponding to each of the extremal states, a* / \/V, 
a*/VV and b#/VV all converge weakly to a complex numbers and 



lim —u) (a?a q ) 

V— »oo V 



lim — uj (a* a ) 



lim —= 

v->oc Jy 



lim —= 



u) (a ) 



u {a Q 



45Q 



4Q(5 + e{q)) 



lim —uj (b*b q ) 

V— »oo V 



lim —= 



u Q (b q 



4(S + e(q))5 



We have again 



p(x) = lim 

V—kx V 



e^ x u (ala p ) 



p + lim — 



fceA* peA* 



e« k -^u (a* k a p ) 



(4.18) 

(4.19) 
(4.20) 

(4.21) 



fceA* peA* p^fc 



The important difference here is that in this model the same boson atoms may condense in 
two states and therefore 



p(x) = p+ lim -2Kc{ 



e iqx uj (a* a q )} 



e^-<^ (a*a 9 )} 



+ lim - V2 3?c{e 

fceA* k^0,q 

+ J im T7 E 23?e ( e ^Kao)}. 



(4.22) 



fceA* fc^0,g 



The last two sums in (|4.22j) vanish in the thermodynamic limit by the same argument as for 
Model 1. However 



lim -u (a*a ) = lim — =u (a*) lim — =u; (a ) = C ^ 0. (4.23) 

V-+oo V * V— »oo yjy l/^oo yjy 

Therefore, the bosons contained in the two condensates may interfere and by virtue of ()4.22j) 
and ()4.23j) this gives the matter-wave grating formed by the quantum interference of the two 
coherent states with different momenta: 

p(x) = p + (Ce iqx + Ce~ iqx ). (4.24) 

Notice that by (J4.23)) and by (|4.24j) there is no matter-wave grating in the Solution 2, when 
one of the condensates (in this case the q- condensate) is empty, see (|2.45|) . 



5 Concluding Remarks 

We conclude this paper with few remarks concerning the importance of the matter-wave 
grating for the amplification of light and matter waves observed in recent experiments. 
It is clear that the absence of the matter-wave grating in Model 1 and its presence in Model 
2 provides a physical distinction between Raman and Rayleigh superradiance. Note first 
that matter-wave amplification differs from light amplification in one important aspect: a 
matter-wave amplifier has to possess a reservoir of atoms. In Models 1 and 2 this is the BE 
condensate. In both models the superradiant scattering transfers atoms from the condensate 
at rest to a recoil mode. 

The gain mechanism for the Raman amplifier is superradiant Raman scattering in a two- 
level atoms, transferring bosons from the condensate into the recoil state pQ 
The Rayleigh amplifier is in a sense even more effective. Since now the atoms in a recoil state 
interfere with the BE condensate at rest, the system exhibits a space matter-wave grating 
and the quantum-mechanical amplitude of transfer into the recoil state is proportional to 
the product of the boson occupation numbers no(n q + 1) for the wave-vectors k = 0,q. 
Each time the momentum imparted by photon scattering is absorbed by the matter-wave 
grating via the coherent transfer of an atom from the condensate into the recoil mode. Thus, 
the variance of the grating grows, since the quantum amplitude for scattered atom to be 
transferred into a recoiled state is increasing [2]-[I], ^2]. At the same time the dressing 
laser beam prepares from the BE condensate a gain medium able to amplify the light. The 
matter-wave grating diffracts the dressing beam into the path of the probe light resulting in 
the amplification of the latter 0. 

In the case of equilibrium BEC superradiance the amplification of the light and the matter 
waves manifests itself in Models 1 and 2 as a mutual enhancement of the BEC and the 
photons condensations, see Solutions 3 in Sections 2 and 3. Note that the corresponding 
formulae for condensation densities for Model 1 (l2~EU|) - (|2~fi2'J) and for Model 2 (|3~TTJ) . (|3~T2J) . 
()3.13|) are identical. The same is true for the boson-photon correlations {entanglements) 
between recoiled bosons and photons, see (|2.64|) . (j3.14|) . as well as between photons and the 
BE condensate at rest: 



r 1 / * , v 1 / *, v Att(\p + e(q) - p)^Xp- p 
v 1 ^ v \ a o- b i/H^, v , P ) = y(aA)Ht* A {w, P ) = ~ 2 , (5-1, 



and for the off- diagonal coherence between recoiled atoms and the condensate at rest: 



r 1/ * x 1 / *\ 4£y(Ap + e(g)-//)(Ap-//) 

As we have shown above, the difference between Models 1 and 2 becomes visible only on 
the level of the wave-grating or spatial modulation of the particle density ([4.24)1 . 
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